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We analyze the spectral function of the single-impurity two-terminal Anderson model at finite 
voltage using the recently developed diagrammatic quantuiu Monte Carlo technique as well as 
perturbation theory. In the (particle-hole-)symmetric case we find an excellent agreement of the 
numerical data with the perturbative results of second order up to interaction strengths U/T ~ 2, 
where T is the transparency of the impurity-electrode interface. The analytical results are obtained 
in form of the nonequilibrium self-energy for which we present explicit formulas in the closed form 
at arbitrary bias voltage. We observe an increase of the spectral density around zero energy brought 
about by the Kondo effect. Our analysis suggests that a finite applied voltage V acts as an effective 
temperature of the system. We conclude that at voltages significantly larger than the equilibrium 
Kondo temperature there is a complete suppression of the Kondo effect and no resonance splitting 
can be observed. We confirm this scenario by comparison of the numerical data with the perturbative 
results. 

PACS numbers: 73.63.Kv, 75.20.Hr, 73.23.-b 



I. INTRODUCTION 

Despite its simplicity the Anderson impurity model 
(AIM)* contains a surprising amount of interesting 
physics. Under equilibrium conditions, when the singly 
occupied dot level lies deep below the chemical poten- 
tial of the fermion continuum and when the Coulomb 
repulsion U prohibits double occupation of the dot, the 
impurity spectral function, i.e. the impurity density of 
states (DOS), is known to develop a sharp Kondo (aka 
Abrikosov-Suhl) resonance which is located at the chem- 
ical potential of the fermion continuum. It is observable 
at relatively low temperatures T < Tk, where the Kondo 
temperature Tk is an estimate for the resonance width. 
This is the essence of the conventional Kondo effect ^^ 

In recent years a quite natural extension of Anderson's 
original idea, namely a setup in which the impurity is 
coupled to two fermion continua, came to attention as 
its experimental realization became feasible^— A partic- 
ularly interesting direction of research is the investigation 
of the nonequilibrium transport as well as of the Kondo 
effect which has been quite successfully approached an- 
alytically, see e . g . Refs. J7H161]. as well as numerically, 
see e. g. Refs. fl7H22l |. Yet another but related research 
direction is concerned with the nonequilibrium transport 
in a pure Kondo model, see e. g. Refs. |23l - [27| . 

In order to induce a finite electric current through a 
quantum dot at least two electrodes with different chem- 
ical potentials are necessary. Contrary to the naive ex- 
pectation that coupling to two fcrmionic continua may 
lead to a two-channel Kondo effect the original resonance 
was predicted to split into two peaks which are positioned 
at the two chemical potentials^ ' 24 ! 25 ' 28 ' 29 Electron trans- 
port between the leads is accompanied by spin-flip pro- 



cesses which essentially break the symmetry necessary to 
drive the system into the two-channel Kondo fixed point. 
On the other hand, exactly those spin-flip processes par- 
ticipate in the RG flow and grow stronger toward low 
energies thereby being responsible for almost perfect ef- 
fective electron transmission (at the Fermi edge in equi- 
librium) across the impurity far below Tk- In the case 
of a finite applied voltage V there is a constant spin- 
flip rate associated with the current flowing through the 
system, so that the Kondo peaks are weaker and broad- 
ened in comparison to the equilibrium situational 3 ^ This 
splitting as well as the broadening of the Kondo peaks 
are difficult to access in the two-terminal setup. There- 
fore a three-terminal approach for the spectral function 
measurement has been put forward ; 31 ! 32 which shortly 
afterwards was implemented experimentally! 33 ! 34 

Although by now a number of studies of multi-terminal 
Kondo /Anderson setups have been conducted, see e. g. 
Refs. 3FJ36| . it is desirable to analyze the problem with 
non- approximative methods in order to cover the full 
Kondo crossover J^- The recently developed diagrammatic 
quantum Monte Carlo (diagMC) approach not only al- 
lows to simulate finite voltage transport but also reli- 
ably works even at zero temperature M-~— In this paper 
the diagMC is applied in order to calculate the impurity 
spectral function of the AIM under nonequilibrium con- 
ditions with special emphasis on the Kondo effect related 
features. The numerical data is then compared with the 
results of the analytical perturbative expansion in inter- 
action strength U. We find that as soon as the bias volt- 
age exceeds the equilibrium Tk the Kondo effect related 
features in the spectral function rapidly deteriorate. By 
comparison of the numerical data with the perturbative 
results we conclude that the large finite voltage has sim- 



(a) 



E 

ED + 



Hi 



u 



eV 




FIG. 1: (a) Level structure of the quantum dot with respect 
to the two fermion continua (leads) at chemical potentials 
Pl/r- The dot level Ed corresponding to single occupation lies 
below Pl/r while the energy level for double occupation Ed + 
U > Pl/r is augmented by U due to Coulomb repulsion, (b) 
Sketch of the setup. For the analytical analysis we consider 
the two terminal setup in which the quantum dot is coupled to 
a left (L) and right (R) lead via electron tunneling oc T L / R = 
r. For the numerical analysis it is advantageous to add a third 
measuring electrode M which is only very weakly coupled to 
the quantum dot and to consider the limit Ta/ — > 0. 



ilar consequences as finite temperature. This explains 
why the Kondo features are so weak and the resonance 
doubling is not observable at all. 

The outline of the paper is as follows: In Sec. lU we 
formulate the problem and identify the most interesting 
energy regimes. The perturbative results for the self- 
energy and impurity DOS, for which we obtain analytical 
expressions in a closed form, are presented. Section ITO1 is 
devoted to the details of the numerical implementation 
of the diagMC scheme. Results are presented in Section 
IV! where the numerical data are compared with the per- 
turbative result and the physical picture is discussed. 



II. FORMULATION OF THE PROBLEM AND 
PERTURBATIVE RESULTS 

We model the system in the canonical way using the 
Hamiltonian 



H — Hq + Hi + H t ■ 



(1) 



Here H contains both the uncoupled impurity level at 
energy e^ which can be subject to a (local) magnetic field 
h, as well as the left/right contacting leads. The latter are 
assumed to be noninteracting fermion continua with field 
operators tp a ,a-, ol = L/R, kept at chemical potentials 
(J-l/r- Here a is the spin index. 



H = ^2 H a [ip a ,o] +^( e D+ HBgcrh/2) d\d a 



(2) 



where d a and d\ are the annihilation/creation operators 
of an electron on the impurity level. As usual, ps is 
Bohr's magneton and g is the Lande factor. Electron 
exchange between the electrodes and the impurity takes 



place locally at x — and is accomplished by 

#t = 7^4V' q , ct (0) + H.c., (3) 



where 7 is the tunneling amplitude between dot level and 
electrode, which we for simplicity assume to be equal 
for both contacts^ Finally, Coulomb repulsion on the 
impurity is taken into account via the last term, 



Hi 



U n-j- nj, , 



(4) 



where n a = d* a d a 



Perhaps the most interesting parame- 
ter regime is the particle-hole symmetric one, when h = 
and cd = —U/2 (aka symmetric Anderson model), for 
which the level structure of the quantum dot is depicted 
in Fig. QJa). Under equilibrium conditions this model is 
solvable by the Bcthc Ansatz^.— This method works per- 
fectly as far as thermodynamic properties are concerned 
but the extraction of single particle quantities, although 
in principle possible, still remains an open issue. The 
most important single particle quantity is the local im- 
purity DOS (spectral function) Pd(w). It is related to the 
Fourier transform of the retarded Green's function (GF) 
of the dot D R (u) via 



where 



Pd {uj) = -2liaD R (u J ). 



D R (t) = -ie(t){{d„(t), 4(o)}) 



(5) 



(6) 



Here {., .} denotes the anticommutator and O(t) is the 
Heavyside step function. In the interacting case, D R (w) 
can be expressed in terms of the self-energy Y, r (oj), 



D R (u) 



1 



tT-E*(w)' 



(7) 



Here T = 2irpo\'f\ 2 is the lead-dot contact transparency 
that depends on the tunneling amplitude 7 and on the 
local DOS in the electrodes po, which we assume to be 
only weakly energy-dependent. For simplicity of nota- 
tion, we will use units for which T = 1 in the following, 
so that r is the unit of energy. 

As has been pointed out in Ref. [43[ the spectral func- 
tion plays an important role especially in nonequilibrium 
transport and the current through the dot is given by the 
Mcir-Wingreen formula^ 



I(V) = G / diL>p d (u>) [n L (io) - Ur(u})] , 



(8) 



where Go = 2e 2 /h is the conductance quantum and 
^L,fl(w) = [exp ((w — pl.r)/T) + 1]~ are the Fermi dis- 
tribution functions in the respective electrodes. First cal- 
culations of the spectral function use perturbation the- 
ory in u£2r— The corresponding series turns out to be 
well-controlled and rapidly converging. Often there is 
no necessity for having the complete energy dependence 





one and + for the anti-timc-ordercd one. 

The lowest order self-energy, c.f. Fig. [2j contains two 
energy integrations and is given by 



FIG. 2: Feynman diagrams for the polarization loop (left), 
and for the self-energy at second order, oc U 2 (right). 



of the spectral function (via e.g. the corresponding self- 
energies) at hand so that an expansion of pd{u) for low 
energies is sufficient. In all other situations there exist 
numerous studies also in nonequilibrium^ i 10 ' 37-41 ' 50 " — 
However, analytical expressions in a closed form for the 
self-energy from the second order onwards do not yet ex- 
ist. 

For technical reasons the calculation of the self-energies 
has a number of advantages. Starting point are the 
impurity GF without interactions. With Hj = the 
Hamiltonian ([T]) is quadratic and can be trivially diag- 
onalizcd even in nonequilibrium^ The noninteracting 

GF (we concentrate on the electron-hole symmetric case) 
read 55 i 56 



D R > A (u>) 

D K {u) 



1 



cj ± i 

2i(n L (uj) +n R (uj) - 1) 
w 2 + l 



(9) 



Just as well we can work with the Keldysh GF D kl (ui) 
which arc related to the GFs in the RAK representation 
([9]) by a simple rotation in Keldysh spaced As usual, the 
Keldysh indices k,l = ± denote the branch of the contour 
on which the final and initial times of the corresponding 
GF in the time domain are taken: — for the time-ordered 



s fcZ 



UJ 



^±D kl ( U -Q)D lk (e)D k \e + Q) 

^D kl ^ + e)U lk (e), (10) 



where 



n' fe (w) = i f — D lk (w + e) D kl (e) 

J 27T 



(11) 



is the generalized (Keldysh) polarization loopi^i An 
explicit calculation of these quantities is tedious but 
straightforward. In the following we give the results only. 
The retarded polarization loop is given by^ 



U r (uj) = H~ (u;)-n-+(u;) 



1 



1 



7T cj 2 + 2iuj 



In 



(1 - ilu) 2 + (T//2) 2 



1 + {V/2f 



(12) 
(13) 



and we find the time-ordered polarisation loop to read 



n-( w ) = n*(| w |) 



-Im 



uj 2 + 2i\uj\ 



In 



\u)\-V/2 + i 
V/2 + i 



e(|7|-H).(i4) 



The remaining Keldysh components can be inferred from 
relation (|12[) and the symmetry properties of II fci , namely 
n+-(w) = II- + (-w) and n+'+(w) = -n~(w)*. With 
these expressions at hand we perform the integral (|TU|> 
and find the retarded self-energy to be given by 



J 



E fl ( W ) = - 




TT 2 L0 1 



(i + uj)(3i + u) 



In 2 



2 x l 2 

n — — In 



i-V/2 
t + V/2 

i-V/2 
i + V/2 

I 



E Li 2 (j9 

a,0=± 



aV/2 - w 
aVj2 + /3i 



^p ( 2i + uj-af3V/2 \ 



(15) 



Here Li2 denotes the dilogarithm function and we have Finally, the spectral function is given via Eqs. ([5j|7|- Its 



introduced 

A(z) = Li 2 (z) - iirsgn[Im(z)} ln(z). (16) 

In equilibrium this expression simplifies considerably, 



__L_^ + 3i Li 2 M 



3/ 



TT 2 (i + u)(3i + uj) 



A(2 -iui) 

4 



.(17) 



behavior as function of energy is depicted in Fig. [S] and 
discussed in section ITvT, 



III. DIAGRAMMATIC MONTE CARLO 
SIMULATION METHOD 

The recently developed refinements of the diagMC 
technique allow to access regimes of arbitrary interac- 



tion strength and therefore make diagMC simulations a 
suitable tool for numerically investigating the AIM with 
large onsite Coulomb repulsion^. Moreover, this ap- 
proach gives access to the transient behavior of the sys- 
tem after sudden switching of the tunneling coupling and 
therefore contributes to the rapidly developing research 
area investiga ting in teraction quenches in quantum dots, 
sec e.g. Rcfs. 27.58 63]. An advantage of diagMC is the 
possibility to explicitly implement the band structure of 
the electrodes. In this regard it can be seen as an ideal 
tool for discussing experimental results. 

In principle, the dot spectral density can be measured 
using its immediate definition ([5]). However, as suggested 
in Refs. |3lll32 |. the spectral function can more efficiently 
be deduced by adding an additional third probe elec- 
trode, which is only weakly coupled to the quantum dot, 
c.f. Fig. [Ub). The conductance of this 'measuring' elec- 
trode (denoted by the subscript 'M' in our notation) is 
given by 



9M{t) = d^ hl{th 



(18) 



where Im and /j,m are the current and chemical potential 
for the third contact, respectively. In the limit of van- 
ishing coupling Tm, the spectral density is related to the 
steady state value (t — > oo) of the conductance b y 31 ' 32 



Pdiv 



M) 



lim 
r M ->o 



r -i 

1 M 



,hm g M (£) 



(19) 



An advantage of this approach is that it corresponds to 
the method used in experimental studies^ 

Following the lines of Refs. J37ll39ll52l ] the current 
through the contact a, with a G {L, M, i?}, in the ab- 
sence of magnetic fields is given by 

I aCT (t) = -4Im^7 CTfc (t)tr{Woa^ fc (t)d CT (t)} , (20) 
feei 

where a aak is the momcntum-A: component of the ^ aa (x) 
field operator and Wo denotes the full initial density op- 
erator of the system (which we assume to factorize, with 
the dot being initially empty). Note that the tunneling 
amplitude jak in Eq. (|20[) is explicitly energy/momentum 
dependent. This is necessary because contrary to the 
analytical calculations of the previous section, where an 
infinite flat band significantly simplifies the calculations, 
the use of a band with a finite width is essential in nu- 
merical approaches for exactly the same reason. An ad- 
ditional speed-up of the diagMC simulations is achieved 
by a smooth switch-on procedure accomplished by a time 
dependent coupling 



lakit) = g a {t)"f ak 



(21) 



where g a (t) interpolates smoothly between at t = and 
1 within the switching time t sw ££(g a (t) = sin 2 [7rf/(2T sw )] 
for t < t sw and g a (t) = 1 for t > r sw ) In the limit of a 
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FIG. 3: Contour-ordered (top) and time-ordered (bottom) 
sequences of tunneling times. 



quasi-continuous distribution of electronic energies in the 
leads, the time-dependent contact transparency reads 

T a (e,t) = 27rp a (e)\ 7ak (t)\ 2 =gl(t)T a (e), (22) 

where /0 Q (e) is the DOS of the leads^ For the numeri- 
cal implementation we choose a flat band with the width 
2e c and temperature-smoothed boundaries^ The corre- 
sponding profile of To^e) is given by 



r«(e) 



\l + ePl*-eJ] \i + e -P(*Uc)] 



(23) 



Time-dependent transport properties can be calcu- 
lated by a diagMC scheme using the conventional 
breakdown of the time evolution operator. Since the 
corresponding formalism has been described in detail 
before , 39 ' 52 we restrict ourselves in the following to pre- 
senting the extensions of the formalism, i.e. the inclusion 
of the probe electrode. 

Expanding the time evolution operator in terms of the 
dot-lead coupling yields a Dyson series where the inte- 
gration variables corresponds to the times at which tun- 
neling between the dot and the lead occurs. The contour- 
ordered kink sequence of the 2n a tunneling times s? of 
spin-cr charges is denoted by 



{s 1 , s 2 , • ■ • , s 2 n } ■ 



(24) 



They reside on the closed real-time contour s : — >• t — >• 
(c.f. Fig. [3}. The corresponding time ordered analogon is 
denoted by 



with 



{S 1 , S 2 , ■ • ■ , S2n„J 1 



for j < n a , 
s°: for j > n a , 



(25) 



(26) 



where n a counts the number of spin-cr tunneling times 
on the forward time axis. Equation (|20[) can then be 



conveniently expressed as 
I a (t) = 



xRe 



with the abbreviation 



2e £ (-l)"t+^ / d? t / d^ 

"4-=° 

{4 a) (4)A(^)^(4.^)} . (27) 



2/ 



ds n 



2/ 



^2n a 



« S 2ti„-1 



ds? (28) 



and the restriction that for a =j" no integration is 
performed over the fixed measurement time s n „+i = 
s rv+ i = t. The influence of the contacts is now sum- 
marized in 



CM) = (-ir^det(5(^)), 
4 a) (? t ) = z^det^^)), 



where 



5y(n) = 



s -+f J- 



(4/-i>4i) f ori<i, 
S + -(4j-i:4) for»>J, 



(29) 



(30) 



:(«)/ 



and SK (sf ) is obtained from Sij(s^) by replacing E fc ' by 
£^ whenever one of the time arguments is equal to the 
measurement time s„ +1 = t. Here, £~ + (E+ - ) denotes 
the dot's lesser (greater) self energy with respect to the 
a contact, and £ fc ' = Y^a^a*^ Using the bandwidth 
profile (f2"B")) we obtain 



g(s)g(s')T a 



X- k ' k (s,s') = 



2/3smh(ir(s-s')//3) 



D kPfj, a 



-ie c (s — s') 



ie c (s— s') 



l_ e -2^ c 



+ 



e -ip, a (s-s') 



-fc^e c 



( e -/s«c _ e -/3M») ( e -/3«= - e^°) 



(31) 



with k = ±. Finally, the contribution from the dot oper- 
ators to Eq. (|2"7|) is given by 



with 



S(«t,S4.) =^ > («t)^ , ( 5 4-) w ( 5 t> 5 4-> n 4-): 



eD(s?-sJ+sf-...)pieD(s?s (T - s 2n„-l+ s 2 



D(s CT ) = e'"^ 1 



and 



(32) 



(33) 



U(s t , s x ) = exp{z[/[rj£ ble (s t , a±) - r^J^, sj.)] } . 

(34) 

Here Tj^ bl denotes the time during which the dot is 
doubly occupied on the forward/backward contour. 



In order to calculate the conductance gu (i) as a func- 
tion of hm hi the limit of vanishing tunneling coupling 
to the third contact (Tm — > 0), we first note that 



-fe.fe 



(s) = -j— -E M ' (s) ex r A/ . (35) 



Next we proceed by collecting all contributions to 
d 



dfj,M 



det(S) det(5 (M) ) 



(36) 



in lowest (linear) order in Tm while neglecting all higher- 
order terms. In the limit Tm — > we find that 



■ dct(S) ex Tm, 



d/i 



dfj,M 

dct(S {M) ) oc T M , 

-det(5 (M) ) = det(S)+0(T 2 M ). 



(37) 
(38) 
(39) 



with 






L^M^lk-nt) foij = n t + l,k<j, 



) dh^ti (4fe-i . *) for j = n t + l,k>j. i 
£ -+ (4k-i> 4j) for 3 ^ n t + !i 3 > k , 



. s+ (4k- 1 . 4j ) for J ^ »t + 1, J > * • 

(40) 
Since [ddet(5)/d//jvf] det(5^ M ') is quadratic in T M while 
det(S)[ddct(S^ M ^) /dfiM] is linear, we arrive at 



dfj,M 



det(5 (M) ) det(S) = det(«§) det(S) + O(r^), 



and the third-terminal conductance is found to read 
0m (4) 



(41) 



lim 

r A f->o r 



,1/ 



-2e 



n^ — 1 



2/ 



(-lft+'H / dst / d$U-l) n + 



2/ 



f det(5(s t )) , , x , - - ,1 , 

xRc | r Af dct ( 5 ^)) ^( ? t.^) | ■ ( 42 ) 



Note that in the zero temperature limit, 

\-k,k f >\ t=o g(js)g(/)Ta i/ia ( s _ s ') 



djl 



2vr 



(43) 



so that gM(t) becomes independent of the bandwidth of 
the third contact. 

A crucial point in the simulation is the attainment of 
the steady state, in which the transport current stops to 
be time-dependent. The further the system is propagated 
in time, starting from the initial time t = when the 
tunneling is switched on, the closer the measured current 
is to the actual steady state value. Unfortunately, this 
time evolution requires a rapidly increasing computation 




Mr 1 ] 



FIG. 4: Spectral function pd{u = Mm) at zero temperature 
and V = 2r for (a) (7 = = en, (im = 2r and (b) U = 4T = 
—en/2, fiM = 4I\ Both panels show data corresponding to 
four different switching times r sw = 0, r _1 , 2I 1- , and 3r _1 
(red, green, blue, and yellow, respectively). The lines are 
guides to the eye only. 
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FIG. 6: Zero temperature dot spectral function for different 
interaction strengths and fixed bias voltage V/T = 2. The 
points are numerical diagMC data and curves are the analyt- 
ical results in the second order in U perturbation theory. The 
slight mismatch of the U — data points and the analytical 
curve is due to finite bandwidth used in the numerical sim- 
ulations (i. e. a bandwidth of 2e c /r = 12 for U/T < 4 and 
2e c /r = 20 for U/T = 8). 



0.3- 



Cl 



0.2- 



0.1- 



u=o, v=o 

U=8T, V=0 
U=8T, V=0.5T 

u=ar, \/=i .sr 
u=8r, i/=2.2r 
U=8r, \/=3r 




FIG. 5: Zero temperature dot spectral function for different 
interaction strengths and bias voltages. The curves shown are 
deduced from the analytical result (|15[) of the perturbative 
calculation. 



time. Still it has been shown in recent work s 37 ' 41 ' 52 that 
in many cases the steady state is indeed accessible with 
moderate numerical effort. As shown in Fig. 21 in our 
case the steady state of gu{t) is typically preceded by 
strongly non-monotonic dynamics. These elongate the 
transient regime and increase the timescale on which the 
steady-state regime is reached. However, this timescale 
can be vastly reduced by adopting a smooth switching 
of the tunneling coupling according to Eq. (p?Tj) , thereby 
extending significantly the parameter regime for which 
the steady state can be reached^ 



IV. RESULTS AND DISCUSSION 

We would like to summarize the results based on the 
second order self-energy first, see Fig. O This approxi- 
mation is known to qualitatively reflect all features one 
expects in the low energy sector. At vanishing voltage 
one obtains the typical trident-shaped curve. The central 
peak in this particular case can be regarded as the pre- 
cursor of the in reality much sharper Kondo (Abrikosov- 
Suhl) resonance 2 - having a width (in equilibrium) which 
can be estimated as 3 . 



V2UT 



' K 



-7rC//(8r) 



(44) 



This formula gives T K - 0.055r for U/T = 8, which 
is much smaller than the actual width of the resonance 
Tk ~ T. The two much wider peaks (shoulders), which 
are located at w w ±U/2 are the Hubbard sidebands. 
Keeping T = and increasing the bias voltage does not 
appear to produce any noticeable qualitative change until 

V hits the threshold of ps T, beyond which the 'Kondo'- 
peak rapidly deteriorates and completely disappears for 

V > r. Interestingly, a very similar destruction of the 
central peak can be observed in equilibrium V = at 
finite temperature T > Tr-. We thus conjecture that 
the effect of the finite voltage might be captured by an 
effective temperature T e g ~ V. 

Further insight is gained through the numerical 
diagMC simulations for the spectral function at finite 
voltage and U. The numerical data, shown in Fig. [51 
turn out to be in an excellent agreement with the analyt- 
ical results using the second order self-energy up to the 



interaction strength U/T = 2. The matching of the per- 
turbative curve and numerical data is still agreeable even 
at U/T = 4. At U/T = 8 one recognizes the very weak 
remnants of the Kondo peak at ui — while the rest of 
the curve has a similar qualitative behavior as the per- 
turbative result and the same order of magnitude. This 
suggests that the complicated collective multi-particle ef- 
fects contained in the full self-energy (i.e. exact in U) 
only have a weak effect and that the main information 
is already accounted for by the lowest order self-energy 
discussed in SecHU This observation is again compatible 
with the existence of an effective temperature. We would 
like to emphasize that T e ff is different from the decoher- 
ence rate discussed, e. g., in Refs. [15lJ25l - 427u30l ] as this is 
not defined in the perturbative regime. 

In the voltage regimes considered, diagMC does not 
produce any evidence for the Kondo peak splitting. The 
relatively high applied voltage V ~ T e g 3> Tk appears 
to induce a widening of the Kondo resonance, which ren- 
ders the observation of the peak splitting impossible. For 
the diagMC approach outlined in Sec. IIII[ smaller volt- 
ages lead to an increase of the timescale over which time- 
dependent transport properties have to be monitored un- 
til a stable stationary regime can be identified. This sig- 
nificantly increases the computational effort necessary to 
extract the steady-state values. Therefore, in the sys- 
tem under consideration and the numerical scheme and 
equipment used, it is not yet possible to give a final an- 
swer to the question of the Kondo peak doubling with 
satisfactorily precision. 

A very interesting issue is the characteristic time scale 
necessary for the Kondo effect to fully develop. In the 



context of a sudd en g ate voltage quench this problem was 
discussed in Ref. |65| . On general grounds one would ex- 
pect that the Kondo peak develops on the time scale 
~ 1/Tk, as Tk is the only energy scale available in the 
system in equilibrium. We observe, however, that even 
in the case of large U — 8T the steady state of the spec- 
tral density at the position of the Kondo peak u> = 
does not take a longer time to establish than for the en- 
ergies outside of the Kondo peak. The only exception 
are the Hubbard subbands uj « ±17/2, where the time 
development appears to be very slow. 

To summarize, we present numerical and analytical re- 
sults for the impurity spectral function of the symmetric 
Anderson model in nonequilibrium at zero temperature. 
We find an excellent agreement of the numerically exact 
diagMC data with the perturbative expansion of the low- 
est order in interaction. For larger U we observe a small 
peak due to the Kondo resonance, which docs not show 
a splitting due to finite applied voltage. 
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